Measurement of Single Electron Event Anisotropy in Au+Au Collisions at y^s/vjv = 

200 GeV 

S.S. Adler,5 S. Afanasiev,i^ C. Aidala,^ N.N. Ajitanand,''^ Y. Akiba,20.38 J. Alexander,^^ Amirikas,!^ 
L. Aphecetche,"'^ S.H. Aronson,'^ R. Averbeck,*'* T.C. Awes,^^ R. Azmoun,'''' V. Babintsev,^^ A. Baldisseri,i° 
K.N. Barish,^ P.D. Barnes,^^ B. Bassalleck,^^ S. Bathe,30 S. Batsouli,^ V. Baublis,^'^ A. Bazilevsky,^^' is 
S. Belikov,i6.i5 Y. Berdnikov,'"' S. Bhagavatula,!^ J.G. Boissevain,^^ H. Borel/° S. Borenstein,^^ M.L. Brooks,^'^ 
D.S. Brown,34 N. Bruner,^^ D. Bucher,30 H. Buesching,30 V. Bumazhnov,^^ G. Bunce,^'^^ J.M. Burward-Hoy,26^44 
S. Butsyk 44 X. Camard 45 J.-S. Chai,!^ p. Chand 4 W.C. Chang,^ S. Chemichenko,!^ C.Y. Chi,^ J. Chiba,20 

M. Chm,9 I.J. Choi,S2 J. Choi,!^ R.K. Choudhury,4 T. Chujo,^ V. Cianciolo,^^ Y. Cobigo,i° B.A. Cole,^ 
P. Constantin,^^ D. d'Enterria,4S G. David, ^ H. Delagrange,45 A. Denisov,^^ A. Deshpande,^^ E.J. Desmond, ^ 
A. Devismes,44 O. Dietzsch,4i O. Drapier,^^ A. Drees,44 K.A. Drees,^ R. du Rietz,^^ A. Dumm,!^ D. Dutta,4 
Y.V. Efremenko,35 K. El Chenawi,49 A. Enokizono,i4 H. En'yo,^^^^^ S. Esuini,48 L. Ewell,^ D.E. Fields,33^39 
F. Fleuret,25 S.L. Fokin,23 B.D. Fox.^s Z. Fraenkel,^! J.E. Frantz,^ A. Franz,^ A.D. Frawley,!^ S.-Y. Fung,^ 
S. Garpman,29'EI T.K. Ghosh,49 A. Glenii,46 G. Gogiberidze,46 M. Gonin,^^ J. Gosset,i° Y. Goto,39 
R. Granier de Cassagnac,^^ N. Grau,^^ S.V. Greene, 4^ M. Grosse Perdekamp,'^^ W. Guryn,*^ H.-A. Gustafsson,^^ 

T. Hachiya,^4 j s. Haggerty,^ H. Hamagaki,* A.G. Hansen, ^"^ E.P. Hartouni,^'^ M. Harvey,^ R. Hayano,^ 
N. Hayashi,38 X. He," M. Heffner,^^ T.K. Hemniick,44 J.M. Heuser,44 M. Hibino,'5° J.G. Hill,!^ W. Holzmann,43 
K. Homma," B. Hong,22 A. Hoover,34 T. Ichihara,38.39 y.V. Ikonmkov,23 K. Imai,24.38 d. Isenhower,i 
M. Ishihara,38 M. Issah,43 A. Isupov," B.V. Jacak,44 W.Y. Jang,^^ Y. Jeong,!^ J. Jia,44 0. Jinnouchi,38 
B.M. Johnson,^ S.C. Johnson,^^ K.S. Joo,3i D. Jouan,36 S. Kametam,^' 5° N. Kamihara,47. 38 j.h. Kang,^^ 
S.S. Kapoor,4 K. Katou,^" S. Kelly,^ B. Khachaturov,^! A. Khanzadeev,37 J. Kikuchi,50 D.H. Kim,3i D.J. Kim,52 
D.W. Kim,i9 E. Kim,42 G.-B. Kim.^s H.J. Kim.^^ Kistenev.^ A. Kivomichi,48 K. Kivovama.32 C. Klein-Boesii " 
H. Kobayashi,38'39 l_ 



ou,^" S. Kelly," B. Khachaturov,^^ A. Klianzadeev,'=" J. Kikuchi,''" D.H. Kim,'^^ D.J. Kim,-^^ 
2 G.-B. Kim,^''^ H.J. Kim,^^ E. Kistenev,^ A. Kiyomichi,48 K. Kiyoyama,32 C. Klein-Boesing,30 
H. Kobayashi,'''''-^'' L. Kochenda,37 V. Kochetkov," D. Koehler,33 T. Kohama," M. Kopytine,44 D. Kotchetkov,^ 
A. Kozlov,5i P.J. Kroon,^ C.H. Kuberg,!'^^ K. Kurita,39 Y. Kuroki,48 M.J. Kweon,^^ Y. Kwon,52 G.S. Kyle,34 
R. Lacey,43 V. Ladygin,!^ J.G. Lajoie,^^ A. Lebedev,^^' ^3 Leckey,44 D.M. Lee,^^ S. Lee,^^ M.J. Leitch,^^ X.H. Li,^ 
H. Lini,42 A. Litvinenko,!^ M.X. Liu,^^ Y. Liu,36 C.F. Maguire,49 Y.I. Makdisi,^ A. Malakhov,!^ V.I. Manko,23 
Y. Mao,^'38 G. Martinez,45 M.D. Marx,44 H. Masui,48 F. Matathias,44 T. Matsumoto,^' P.L. McGaughey,^^ 
E. Melnikov,i5 F. Messer,44 Y. Miake,48 J. Milan,43 T.E. Miller,49 A. Milov,44^5i S. Mioduszewski,^ R.E. Mischke,^^ 

G. C. Mishra,i3 J.T. Mitchell,^ A.K. Mohanty,4 D.P. MoiTison,'^ J.M. Moss,^^ F. Muhlbacher,44 D. Mukhopadhyay,^^ 

M. Muniruzzaman,^ J. Murata,38'39 Nagamiya,^^ J.L. Nagle,^ T. Nakamura,i4 g k. Nandi,'' M. Nara,48 
J. Newby,46 P. Nilsson,29 A.S. Nyanin,23 J. Nystrand,^^ E. O'Brien,^ C.A. Ogilvie,!^ H. Ohnishi,5'38 i_d. Ojha,49.3 
K. Okada,38 M. Ono,48 V. Onuchin,!^ A. Oskarsson,^^ I. Otterlund,^^ K. Oyama,* K. Ozawa,^ D. Pal,^! 
A.P.T. Palounek,27 V. Pantuev,44 V. Papavassiliou,34 J. Park,42 A. Parmar,33 S.F. Pate,34 T. Peitzniann,30 
J.-C. Peng,2^ V. Peresedov,!"^ C. Pinkenburg,^ R.P. Pisani,^ F. Plasil,35 M.L. Purschke,^ A.K. Purwar,44 
J. Rak,i6 1. Ravinovich,5i K.F. Read,35.46 m. Reuter,44 K. Reygers,30 V. Riabov,37^4o y. Riabov,37 G. Roche,^^ 
A. Romana,25 M. Rosati,!^ p. Rosnet,^^ S.S. Ryu,^^ M.E. Sadler,^ N. Saito,38'39 T. Sakaguchi,^' 5° M. Sakai,32 
S. Sakai,48 V. Samsonov,37 L. Sanfratello,33 R. Santo,30 H.D. Sato,24^38 Sato,5^48 gawada,20 Y. Schutz,45 
V. Semenov,i5 R. Seto,^^ M.R. Shaw,i'27 T.K. Shea,^ T.-A. Shibata,47^ 38 Shigaki,i4.20 x. Shiina,27 C.L. Silva,4i 
D. Silvermyr,27^29 g^ ^^^22 Q p gjngii^s y. Singh,3 M. Sivertz,^ A. Soldatov,!^ R.A. Soltz,26 W.E. Sondheim,27 
S.P. Sorensen,46 I.V. Sourikova,^ F. Staley,i° P.W. Stankus,35 E. Stenlund,29 M. Stepanov,34 A. Ster,2i 
S.P. Stoll,^ T. Sugitate," J.P. Sullivan,27 E.M. Takagui,4i A. Taketani,38.39 M. Tamai,^'' K.H. Tanaka,20 
Y. Tanaka,32 K. Tamda,38 M.J. Tannenbaum,^ P. Tarjan," J.D. Tepe,i'27 T.L. Thonias,33 J. Tojo,24.38 

H. Torii,24.38 j^ g^ Towell,^ I. Tserruya,^! H. Tsuruoka,48 S.K. Tuli,3 H. Tydesjo,^^ N. Tyurin,!^ H.W. van Hecke,^^ 
J. Vclkovska,5'44 M. Velkovsky,44 V. Veszpremi," L. Villatte,46 A. A. Vinogradov,23 M.A. Volkov,23 E. Vznuzdaev,37 

X.R. Wang,i3 Y. Watanabe,38' 39 g Wliite,^ F.K. Wolin,!^ C.L. Woody,^ W. Xie,^ Y. Yang,^ A. Yanovich,!^ 
S. Yokkaichi,38.39 G,R, Young,35 I.E. Yushmanov,23 W.A. Zajc,3'[l| C. Zhang,^ S. Zhou,^ S.J. Zhou,^! and L. Zolin" 

(PHENIX Collaboration) 

^Abilene Christian University, Abilene, TX 79699, USA 
^Institute of Physics, Academia Sinica, Taipei 11529, Taiwan 
^Department of Physics, Banaras Hindu University, Varanasi 221005, India 
'^Bhabha Atomic Research Centre, Bombay 400 085, India 
Brookhaven National Laboratory, Upton, NY 11973-5000, USA 
^University of California - Riverside, Riverside, CA 92521, USA 



2 



China Institute of Atomic Energy (CIAE), Beijing, People's Republic of China 
^ Center for Nuclear Study, Graduate School of Science, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 
'•'Columbia University, New York, NY 10027 and Nevis Laboratories, Irvington, NY 10533, USA 
'"Dapnia, CEA Saclay, F-91191, Gif-sur-Yvette, France 
^'Debrecen University, H-4010 Debrecen, Egyetem ter 1, Hungary 
^'''Florida State University, Tallahassee, FL 32306, USA 
"Georgia State University, Atlanta, GA 30303, USA 
''^Hiroshima University, Kagamiyama, Higashi- Hiroshima 739-8526, Japan 
Institute for High Energy Physics (IHEP), Protvino, Russia 
^'^lowa State University, Ames, I A 50011, USA 
^'^ Joint Institute for Nuclear Research, 14-1980 Dubna, Moscow Region, Russia 
KAERI, Cyclotron Application Laboratory, Seoul, South Korea 
^'^ Kangnung National University, Kangnung 210-702, South Korea 
^'^KEK, High Energy Accelerator Research Organization, Tsukuba-shi, Ibaraki-ken 305-0801, Japan 
KFKI Research Institute for Particle and Nuclear Physics (RMKI), H-1525 Budapest 114, POBox 49, Hungary 

^^Korea University, Seoul, 136-701, Korea 
'^''Russian Research Center "Kurchatov Institute", Moscow, Russia 
^''Kyoto University, Kyoto 606-8502, Japan 
^'^ Laboratoire Leprince-Ringuet, Ecole Poly technique, CNRS-IN2P3, Route de Saclay, F-91128, Palaiseau, France 
''^Lawrence Livermore National Laboratory, Livermore, CA 94550, USA 
'''Los Alamos National Laboratory, Los Alamos, NM 87545, USA 
^^LPC, Universite Blaise Pascal, CNRS-IN2P3, Clermont-Fd, 63177 Aubiere Cedex, France 
'"^''Department of Physics, Lund University, Box 118, SE-221 00 Lund, Sweden 
'^"institut fiir Kemphysik, University of Muenster, D-48149 Muenster, Germany 

Myongji University, Yongin, Kyonggido 449-728, Korea 
'''^Nagasaki Institute of Applied Science, Nagasaki-shi, Nagasaki 851-0193, Japan 
'^^ University of New Mexico, Albuquerque, NM 87131, USA 
^^New Mexico State University, Las Cruces, NM 88003, USA 
'^^Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA 
^^IPN-Orsay, Universite Paris Sud, CNRS-IN2P3, BPl, F-91406, Orsay, France 
''PNPI, Petersburg Nuclear Physics Institute, Gatchina, Russia 
^^RIKEN (The Institute of Physical and Chemical Research), Wako, Saitama 351-0198, JAPAN 
^^RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973-5000, USA 
^"St. Petersburg State Technical University, St. Petersburg, Russia 
Universidade de Sao Paulo, Instituto de Fisica, Caixa Postal 66318, Sao Paulo CEP05315-970, Brazil 
''^System Electronics Laboratory, Seoul National University, Seoul, South Korea 
^''Chemistry Department, Stony Brook University, SUNY, Stony Brook, NY 11794-3400, USA 
'^'^ Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, NY 11794, USA 
^^SUBATECH (Ecole des Mines de Nantes, CNRS-IN2P3, Universite de Nantes) BP 20722 - 44307, Nantes, France 

University of Tennessee, Knoxville, TN 37996, USA 
'^'^ Department of Physics, Tokyo Institute of Technology, Tokyo, 152-8551, Japan 
Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
Vanderbilt University, Nashville, TN 37235, USA 
"'"Waseda University, Advanced Research Institute for Science and 
Engineering, 17 Kikui-cho, Shinjuku-ku, Tokyo 162-0044, Japan 
Weizmann Institute, Rehovot 76100, Israel 
^^Yonsei University, IPAP, Seoul 120-749, Korea 
(Dated: February 8, 2008) 

The transverse momentum dependence of the azimuthal anisotropy parameter V2, the second 
harmonic of the azimuthal distribution, for electrons at mid-rapidity {\rj\ < 0.35) has been measured 
with the PHENIX detector in Au+Au collisions at ^snn — 200 GeV. The measurement was made 
with respect to the reaction plane defined at high rapidities {\ri\ — 3.1 — 3.9). From the result we 
have measured the V2 of electrons from heavy flavor decay after subtraction of the V2 of electrons 
from other sources such as photon conversions and Dalitz decay from light neutral mesons. We 
observe a non-zero single electron V2 with a 90 % confidence level in the intermediate pr region. 
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I. INTRODUCTION 

The azimuthal anisotropy of particle emission is a pow- 
erful tool to study the early stage of ultra-relativistic 
nuclear collisions. The spatial anisotropy in the initial 
stage of non-central nucleus-nucleus collisions is trans- 
ferred into momentum anisotropy in the final state. The 
azimuthal anisotropy is defined by 

^ =iVo|l + X^2w„cos(n((^-*j^.P.))| , (1) 

where A'o is a normalization constant, (j) is the azimuthal 
angle of particles, and '^r.p. is the direction of the nuclear 
impact parameter ("reaction plane") in a given collision. 
The harmonic coefficients, w„, indicate the strength of the 
n*'' anisotropy. The azimuthal anisotropy parameter V2 
(the second harmonic coefficient of the Fourier expansion 
of the azimuthal distribution) may be especially sensitive 
to the early pressure . The transverse momentum (j>t) 
dependence of has been measured for identified par- 
ticles at RHIC 0, IE IE • Previous measurements are 
limited to hadrons made of light quarks. These results 
show a clear mass dependence of V2, which is well repro- 
duced by a hydrodynamical calculation Q in the low px 
region {px < 2 GeV/c). The agreement is considered as 
evidence that the collective motion develops in the very 
early stages of the reaction. It is also observed that V2 as 
a function of px scales via the coalescence prescription, 
that is, V2/n as a function of px/n is universal, where n 
is the number of valence quarks plus valence anti-quarks. 
This scaling behavior is consistent with the prediction of 
the quark coalescence model, which assumes a finite V2 of 
quarks |^ . This suggests that the V2 already develops in 
the partonic phase for hadrons made of light quarks. In 
addition, if the V2 of heavy quarks is non-zero, it would 
support partonic level thcrmalization and very high den- 
sity at the early stage of the collisions. 

Electrons are a useful tool to study the production of 
heavy quarks such as charm quarks. In the PHENIX ex- 
periment at RHIC, transverse momentum spectra of sin- 
gle electrons have been measured in Au-|-Au collisions at 
y^s^N = 130 GeV and 200 GeV The results are 
consistent with that expected from semileptonic charm 
decays in addition to decays of light mesons and photon 
conversions On the other hand, electrons originating 
from semileptonic decays of D mesons have a significant 
angular deviation from the original D meson direction. 
The effect for V2 has been studied in and |22|. The 
results suggest that the effect is not significant for the 
decay electron V2, and the electron V2 reflects the V2 of 
D meson. Therefore the single electron V2 measurement 
is a useful method for studying open charm V2- 

Currently the single electron spectra from PHENIX 
are consistent with two opposing scenarios: (1) initial 
perturbative QCD charm production without final state 
interactions and (2) complete thermal equilibrium for 



charmed hadrons jlOj. Therefore the measurement of 
the azimuthal anisotropy of electrons from semileptonic 
charm decays could give us important new information 
regarding the charm dynamics in high-energy heavy-ion 
collisions. The measurement is also important for the 
quark coalescence model due to the large difference be- 
tween the charm quark and light quark masses, and for 
the prediction of V2 for the J/tp and the D meson, which 
contain charm quarks. 

In this paper, we present the first measurement of the 
single electron V2, which is expected to reflect the heavy 
flavor azimuthal anisotropy, below 4 GeV/c with respect 
to the reaction plane in Au-|-Au collisions at ^/sJfN = 
200 GeV. The single electron V2 was measured by sub- 
tracting from the inclusive electron V2 the V2 of electrons 
from other sources such as photon conversions and Dalitz 
decays from light neutral mesons. 
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FIG. 1; PHENIX experiment configuration in Run2. Top: 
Cross section perpendicular to the beam pipe. Bottom: East 
side view of the cross section along the beam pipe. 



II. DATA ANALYSIS 

About 16M minimum bias events in RHIC-Run2 
(2001) for ^/sNN = 200 GeV are used in this analysis 
after a vertex cut is applied (|zvcrtcx| < 20 cm). In this 
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FIG. 2: The configuration of the 64 PMTs of each BBC is 
shown. From the hit position of the particles in each PMT, 
the azimuthal angle is calculated. 



section we present a brief overview of the PHENIX detec- 
tors ^3 used in this analysis and then present details of 
event selection, electron identification and reaction plane 
determination. 



A. Overview of PHENIX detector 

PHENIX consists of four spectrometer arms (central 
arms and muon arms) and a set of global detectors. The 
central arms are located east and west of the interaction 
region at mid-rapidity. The muon arms are located to 
the north and south at forward rapidity. Fig. ^ shows 
the configuration of the central arms in Run2. 

The global detectors consist of the Beam-Beam- 
Counters (BBCs) and the Zero-Degree Calorimeters 
(ZDCs). These detectors provide the time of the Au-f Au 
collision, the collision vertex, the event trigger and the 
collision centrality. In this analysis the BBCs are also 
used to determine the reaction plane. The BBCs are 
installed on the North and South sides of the collision 
point along the beam axis. Each BBC is placed 144 cm 
from the center of the interaction region and surrounds 
the beam pipe. This corresponds to a pseudorapidity 
range from 3.1 to 3.9 over the full azimuth. Each BBC is 
composed of 64 elements; the configuration is shown in 
Fig. 121 A single BBC element consists of a one-inch diam- 
eter mesh dynode photomultiplier tube mounted on a 3 
cm long quartz radiator. The ZDC is a hadron calorime- 
ter and measures the energy of spectator neutrons. The 
ZDCs are located 18 m downstream and upstream along 
the beam axis, and each ZDC covers 2 mrad of forward 
angular cone, corresponding to rj >6.0. 

The central arms are designed to track particles emit- 
ted from collisions, identify charged particles and recon- 
struct invariant masses. The central arms each cover the 
pseudorapidity range \r]\ < 0.35 and 90° in azimuthal an- 
gle. The central arms consist of several subsystems. In 
this analysis drift chambers (DCs), pad chambers (PCs), 



ring imaging Cherenkov counters (RICHs) and electro- 
magnetic calorimeters (EMCals) are used. The DCs are 
located between 2.0 and 2.4 m from the beam axis on 
each central arm and measure charged particle trajecto- 
ries in the r-0 plane. The central arms have three layers 
of PCs, which are multi-wire proportional chambers. The 
PCs are located at 2.4 m (PCI), 4.2 m (PC2) and 5.0 m 
(PC3) from the beam axis. PCI and PC3 are installed 
in each central arm, but PC2 is installed only in the 
west arm. The PC measures 3-D space points along the 
straight line particle trajectories. A RICH, the primary 
detector for electron identification, is installed in each 
central arm. The RICH consists of a gas vessel, a thin 
reflector and a photon detector consisting of an array of 
PMTs. During Run2 CO2 was used as the Cherenkov ra- 
diator so only pions with p > 4.7 GeV/c emit Cherenkov 
light in the RICH. The EMCal is used to measure the 
spatial position and energy of electrons and photons. It 
covers the full central arm acceptance of 70° < 9 < 110° 
with each of the two walls subtending 90° in azimuth. 
One wall is comprised of four sectors of Pb-scintillator 
sampling calorimeter, and the other has two sectors of 
Pb-scintillator and two of Pb-glass Cherenkov calorime- 
ter. The Central Magnet (CM) provides a magnetic field 
around the interaction vertex that is parallel to the beam. 
The CM allows momentum analysis of charged particles 
in polar angle range from 70° < 9 < 110° and provides a 
field integral of about 0.8 Tesla-meters p^. 



B. Event selection 

The event selection was done with the BBC and the 
ZDC in this analysis. The minimum bias trigger requires 
a coincidence between north and south BBC signals. The 
trigger included 92.2+^-^n% of the 6.9 barn Au-|-Au in- 
elastic cross section [lj|. The event centrality is deter- 
mined by combining information on spectator neutrons 
measured by the ZDC and the charge sum information 
measured by the BBC. The collision vertex point along 
the beam line is determined by the timing difference of 
the two BBCs. We required IzvortcxI < 20 cm for this 
analysis. 



C. Cliarged particle selection and electron 
identification 

Charged particle tracks are reconstructed by the DC 
and the first pad chamber plane (PCI) installed in each 
central arm together with the collision vertex determined 
by the BBC [lj|- In order for a reconstructed track to 
be selected, the track projection to the EMCal and the 
position of the associated hit in the EMCal must match 
within 2 standard deviations. The electron candidates 
are required to have at least three associated hits in the 
RICH that pass a ring shape cut and are also required to 
pass a timing cut. To reduce background from hadrons 
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FIG. 3: {E - p)/p/a distribution. We require -2a < [E - 
p)/p < 3(T to reduce background from hadrons and piioton 
conversions far from the vertex. A background of less than 
10 %, caused by accidental association of tracks with RICH 
hits, remains. 



and photon conversions far from the vertex, energy is 
measured in the EMCal, and momentum matching {E/p) 
is required. Electrons deposit all of their energy in the 
EMCal; therefore the E/p is approximately 1.0. In this 
analysis we require —2a < {E — p)/p < 3cr to reduce 
background. Fig. |2| shows the {E — p)/p/a distribution. 
Here the a means a standard deviation of {E — p) /p. A 
background of less than 10 % remains, caused by acci- 
dental association of tracks with RICH hits. The back- 
ground level is estimated by an event mixing method and 
is subtracted when we calculate the electron V2- 



D. Reaction plane determination 

In this analysis the values of V2 are calculated by the 
reaction plane method, which measures the azimuthal 
angle of the particle emission with respect to the reaction 
plane . The azimuthal angle of the reaction plane for 
the n*'* harmonic is determined by |l6j | 



~ tan 



_i Y^jWi sin(n0t 
Y,i Wi cos{n(f>i 



/n, 



(2) 



where is the azimuthal angle of each particle used 
in the reaction plane determination and Wi is the corre- 
sponding weight. The azimuthal angle distribution of the 
particle emission measured with respect to the reaction 
plane can be written as Eq.^ Due to finite reaction plane 
resolution, coefficients in the Fourier expansion of the 
azimuthal distribution with respect to the "measured" 
reaction plane {v^'^'^') are smaller than coefficients mea- 
sured with respect to the "real" reaction plane (u„). The 
resolution correction necessary for v^°^^ is given by; 



(3) 



where u„ is the real coefficient and CTi,,, is the reaction 
plane resolution for the rt*'' harmonic. The cr^„ is es- 
timated with data using a formula shown in Ref. |l6j| . 
The value of v'^'^^^- is obtained by fitting the azimuthal 
distribution (relative to the reaction plane) with 



— =iVo(l-f 2«r"-cos(n0)), 



(4) 



where Nq and v^'^'^- are fitting parameters. We can also 
calculate v^'^'^^- directly by 



= (cos(n0)). 



(5) 



In this analysis the V2 is estimated by using the re- 
action plane found from the second harmonic (n = 2), 
since better accuracy of w„ is obtained by using the same 
harmonic's reaction plane T^. The reaction planes are 
determined by using both BBCs. In the PHENIX exper- 
iment the reaction plane is also determined by using the 
central arm detectors. One of the key issues of the reac- 
tion plane determination is non-flow effects such as jets, 
resonance decays and HBT. Since each BBC is roughly 
three units of pseudorapidity away from the central arms, 
it is expected that the non-flow effects are smaller there 
than in the central arm detectors 4]. 

Using the BBC information the reaction plane is mea- 
sured by 



1 g»sin(20,) 



tan 



/2 



(6) 



where 4>i is the azimuthal angle of each PMT and qi is 
the charge information of each PMT. Due to the ran- 
dom distribution of the impact parameter direction in 
collisions, the reaction plane should have an isotropic 
azimuthal distribution. Because of the possible az- 
imuthal asymmetries in the BBC response, however, 
the measured reaction plane distribution is anisotropic. 
In this paper, we use the following two step method 
to correct the reaction plane. First the distribution 
of ^^^j^ sin(20i) and ^^^j^ cos(20i) are recentered 
by subtracting (y^^^^ sin(20i)) and (X^^^i cos(20i)) 
over all events jl6j: 



^^^-1 <li sin(20j) - Qi sin(2(/),)) 'j 



* C0S(2(/|,;) - C0S(2(/|,)) 



(7) 

This method does not remove higher harmonic compo- 
nents of the determined reaction plane, so we apply an 
additional correction method this method flat- 

tening the reaction plane is accomplished by using a shift 
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fiat 



nipo 



bs. 



(8) 



where Aijj is the correction factor for the reaction plane. 
Atp is determined by 



' An cos(2nV'obs.) + -S« sin(2nV'obs.)- (9) 
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An and i?„ are defined by requiring the n*^ Fourier mo- 
ment of the new reaction plane (-0**^*) to vanish. 



2 

An = (sin(2n'(/'obs.)) 

n 

2 

Bn = -(cos(2n'(/'obs.))- 
n 



(10) 
(11) 



Since the reaction plane depends on collision centrality 
and z vertex, the reaction planes are divided into 40 sam- 
ples (20 centrality bins and 2 vertex bins), and these cor- 
rections are determined independently for each sample. 
A combined reaction plane, which is defined by weighted 
averaging the reaction plane angles obtained by the south 
side BBC and the north side BBC, is used to measure the 
V2 in this analysis. The resolution of the combined reac- 
tion plane is estimated by using Eq. 11 in [l6]. Figure 01 
shows the centrality dependence of the resolution (cuj). 
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FIG. 4: Centrality dependence of the combined reaction 
plane resolution that is determined by the BBC. 



the background are shown in the right panel, and the 
yields are shown as small symbols in the left panel. The 
electron V2 are measured after subtraction of this back- 
ground {dN^'""''' / d(j)) from electrons that are identified 
by the RICH {dN^"-"" /d<j)): 



dN" 



dN" 



dN""-"- 



(12) 



The transverse momentum dependence of V2 for elec- 
trons after subtracting background is shown in Fig. 
The statistical errors are shown as vertical lines in the 
figure. The 1 a systematic uncertainties are shown as 
vertical bands. The systematic uncertainties include the 
systematic uncertainty of the reaction plane determina- 
tion and electron identification. The systematic uncer- 
tainty of the reaction plane determination is about 5%. 
The uncertainty was estimated by measuring V2 with re- 
action plane which was determined by the North side 
BBC, the South side BBC, and a combination of the 
North and South sides. The systematic uncertainty from 
electron identification was estimated by measuring elec- 
tron V2 with several different sets of electron identifica- 
tion cuts. A comparison with V2 for charged pion Q is 
also shown in Fig. At low pr {pt < 1-0 GeV/c), the 
electron V2 is larger than the V2 of the pion. In this re- 
gion electrons come mainly from tt*^ decays, directly from 
the Dalitz decays or indirectly from photon conversions. 
Due to the fact that the decay angle of the 7r° decay is 
small, the electron has about the same azimuthal angle 
as the parent tt'^, while at the same time the electron px 
is smaller than the tt*' pr- Therefore the electron V2 at a 
given Pt corresponds to the larger V2 of the tt^ at higher 
Pt- The V2 of charged pions is consistent with that of 
neutral pions [T9| . therefore the inclusive electron V2 is 
higher than the tt V2- 



E. Inclusive electron V2 
III. RESULTS AND DISCUSSION 

In this section we present a method to calculate the 
single electron V2 from inclusive electrons and show its 
transverse momentum dependence. 

The azimuthal distributions of electrons relative to the 
reaction plane are shown in the middle panel of Fig. \S[ 
The distributions are overlaid by shifting them on the 
vertical axis such that the spacing between each is equal. 
Each symbol represents the measured pt region indicated 
in the left panel, which shows the raw yields of the each 
distribution with large symbols. As described in section 
III CI less than 10% background remains from acciden- 
tal RICH associations. The azimuthal distributions of 



A. Heavy flavor electron 112 

The inclusive electron sample has two components: 
(1) "non-photonic" - primarily semi-leptonic decays of 
mesons containing heavy (charm and bottom) quarks, 
and (2) "photonic" - Dalitz decays of light neutral mesons 
(ttq, 77, 77', uj and (p) and photon conversions in the detec- 
tor material . The azimuthal distribution of electrons 
(dN^/dcj)) is the sum of the azimuthal distributions of 
photonic electrons {dN'^ /dcj)) and non-photonic electrons 
(dA^"°"-^/d0): 

dNe _ dN2 dN"^""^-^ 
14 ~ ~d4^ d4 ■ 

The second harmonic of the Fourier expansion of each 
azimuthal distribution is defined according to 
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FIG. 5: Azimuthal distributions relative to the reaction plane of electrons (middle panel) and the background (right panel) 
in each momentum bin are shown. The distributions are overlaid by shifting them on the vertical axis such that the spacing 
between each is equal. The raw yields of each distribution are shown the left panel. 



Ne (1 + 2v2, cos(20)) ^ N2 (1 + 2vl^ cos(20)) + Ar»°""7 (i + 2v'!^°'"~^ cos(20)) 

= (iV7 + 1 + 2 - \,^'^non-"' C0s(2<^) , (14) 



where is the V2 of inclusive electron, is the V2 of 
the photonic electrons and V2°^~'^ is the V2 of the non- 
photonic electrons. From Eq. ^1 the inclusive electron 
f 2 is given by: 



V2, 



N2vl + N^on-^^non-^ 
— + (1 — r)V2 , 



(15) 



where r is defined as r = 1/(1 + Rnp). Rnp is the ra- 
tio of the number of non-photonic electrons to photonic 



electrons (7V"°"~'''/A'7). We experimentally determined 
the ratio from analysis of special runs in which addi- 
tional photon converter was installed. The details of the 
method are described in ,21], and the measured ratio is 
shown in Fig. [7| The increase in the number of non- 
photonic electrons is consistent with that expected from 
semileptonic charm decays From Eq.^]?;™""''^ can 
be expressed as 



V2, 



1 



(16) 



The dominant sources of photonic electrons are photon 
conversions and Dalitz decays from tt" I^. In addition, 
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we also took into account electrons from r/ decays when 
calculating photonic electron V2- We assumed that the 
contributions from rj decays is 17 % taken into account 
^/■""o = 0.45 [231 . The other sources are ignored when 
calculating the photonic electron V2 due to their small 
contribution. The decay electron V2 from decay electrons 
of 7r° and 77 were calculated by Monte Carlo simulation. 
The transverse momentum dependence of the tt" V2 was 
obtained from the measured ttO v2 M {pt > 1.0 GeV/c) 
and the measured charged tt V2 (pt < 1.0 GeV/c). Both 
measurements were used since the tt^ V2 has been mea- 
sured only above 1.0 GeV/c, and both V2 measurements 
are consistent at intermediate px (1.0 < pr < 3.0). The 
measured spectra |l8j were used to give the input 
transverse momentum spectrum. The transverse momen- 
tum dependence of the rj V2 was taken to be the same as 
for the kaon V2. The transverse momentum spectrum of 



rj was approximated by assuming rriT scaling of tt'' spec- 
tra. The photonic electron V2 that is calculated from the 
results is shown as the dashed line in the left panel of 
Fig. |S1 The middle dashed line is the mean value of the 
photonic electron V2 and the upper and lower dashed lines 
show the 1 (J systematic uncertainty. The systematic un- 
certainty of the photonic electron V2 was estimated from 
the statistical error and the systematic error of the mea- 
sured parent V2. If the non-photonic electron V2 is zero, 
that means the V2 of the parent particle, such as a I? 
meson, is zero. Additionally, the inclusive electron V2 is 
the same as that of the scaled photonic electron (rw2^) 
from Eq. ^] The scaled photonic electron V2 is shown 
as the solid line in the left panel in Fig. |H1 At interme- 
diate Pt (1.0 < Pt < 1-5) the electron V2 is higher than 
rv2'' ■ This might suggest that the non-photonic electron 
has non-zero V2 at intermediate pt- The details of this 
discussion are presented in the next section. 

Background from kaon decays {K —fTrev) remain in 
the non-photonic yield. The contribution of kaon decays 
to the non-photonic yield, shown in Fig.Qas dashed line, 
is 18 % ai Pt — 0.4 GeV/c and decreases rapidly to less 
than 6 % for = 1 GeV/c |23|. The transverse momen- 
tum dependence of the kaon V2 has been measured up to 
3.0 GeV/c and that of the V2 has been measured up 
to 6.0 GeV/c The kaon and V2 are consistent up 
to 3.0 GeV/c, and the quark coalescence model predicts 
that these two meson V2 values are the same. Therefore, 
kaon and K'g V2 were combined as input for the kaon 
V2. The transverse momentum spectrum of kaons was 
obtained from measured kaon spectra up to 2.0 GeV/c. 
In the high pT region we used scaled tt" spectra and as- 
sumed that the shapes of the kaon spectra were the same 
as the TT*^ spectra, which arc matched with measured kaon 
spectra around 2.0 GeV/c. 

The main source of the non-photonic electrons is 
semileptonic decays of heavy flavor (charm and beauty) . 
Therefore the non-photonic electron V2 that was obtained 
by subtracting photonic electron and kaon decays from 
inclusive electrons should be heavy flavor electron V2 , 
which reflects the azimuthal anisotropy of heavy quarks. 
The result of the heavy flavor electron V2 is shown in the 
right panel of Fig. [S] The vertical lines are the statistical 
errors that are propagated from the statistical errors of 
the inclusive electron V2 shown in Fig. |^ The la sys- 
tematic uncertainty of heavy flavor electron V2 is shown 
as bands. The systematic uncertainty includes the sys- 
tematic uncertainty of the reaction plane, the measured 
inclusive electron V2 (w/o R.P.), the photonic electron 
V2 (w/o R.P.) and Rnp- The systematic uncertainty of 
Rnp is the quadratic sum of the statistical and system- 
atic errors because Rnp is measured with different data 
set. There are two categories of uncertainty: Type A 
is a point-to-point error uncorrected between pT bins, 
and type B is a common displacement of all points by 
the same factor independent of pt- The total systematic 
uncertainty is calculated by propagating the errors on 
the individual quantities that enter into Eq. ^| Tabled 
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TABLE I: The relative systematic uncertainty of heavy quark electron 112. All errors are given in percent. 



Pt range 


Sys. error bound 


inclusive e V2 (%) 


photonic e V2 (%) 


Rnp (%) 


R.P. (%) 


Total (%) 


0.4 <PT < 1.0 


lower 


< 32 


< 26 


< 21 


4.5 


< 42 




upper 


< 63 


< 21 


< 21 


4.5 


< 70 


1.0 <PT < 1.75 


lower 


< 25 


< 21 


< 14 


4.5 


< 36 




upper 


< 67 


< 17 


< 15 


4.5 


< 70 


1.75 <pt < 4.0 


lower 


280 


78 


190 


4.5 


340 




upper 


220 


64 


280 


4.5 


360 


Type 




A 


A 


A 


B 





shows the relative systematic uncertainty of heavy quark 
electron V2- 

From the result we calculated the confidence level for 
a non-zero V2- We assumed that the data of measured 
heavy flavor electron V2 follows a Gaussian distribution 
and the a was obtained by calculating quadratic sum of 
the statistical and systematic errors of the heavy flavor 



electron V2 assumed these errors are independent. In the 
intermediate Pt region (1.0 GeV/c < px < 1.75 GeV/c), 
the confidence level is 90 %, suggesting that the measured 
heavy flavor electron V2 has a non-zero V2- 

Assuming the quark coalescence mo del, decay electron 
V2 from D mesons has been predicted In the model 
D mesons are formed from charm quark coalescence with 
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FIG. 9: Comparison of the heavy flavor electron V2 with two different charm flow scenarios from [TH . The sohd hne corresponds 
to no rescattering of the initially produced charm quarks (without flow), while the dashed line reflects the effect of complete 
thermalization (with flow). 



thermal light quarks at hadronization. For charm quark 
momentum spectra, two extreme scenarios are consid- 
ered. The first scenario assumes no reinteractions after 
the production of charm-anticharm quark pairs in initial 
state hard processes (calculated from PYTHIA). The sec- 
ond scenario assumes complete thermalization with the 
transverse flow of the bulk matter. FigureElshows a com- 
parison of the heavy flavor electron V2 with decay elec- 
trons from D mesons in the "no reinteraction" scenario 
as a solid line, while the dashed line reflects the "thermal- 
ization" scenario. Due to large systematic and statistical 
uncertainty of the current measurement, neither scenario 
is excluded by this single electron V2 measurement. 



IV. SUMMARY 

In summary, we have measured the elliptic flow, V2^ 
of single electrons from heavy flavor decay. This single 
electron V2 is produced by subtracting the V2 of electron 
sources such as photon conversion from the V2 of inclu- 
sive electrons measured with the PHENIX detector in 
Au-|-Au collisions at ^/snn — 200 GeV with respect to 
the reaction plane defined at high rapidities (|ry| = 3 — 4). 
The measured heavy flavor electron V2 is nonzero with a 
90 % confidence level. Two model calculations from 
assume extremely different scenarios: either no reinterac- 
tion of the initially produced charm quarks or complete 
thermalization with the bulk matter. Both of these cal- 
culations are consistent within errors with the measured 
heavy flavor electron V2. 



High luminosity Au-f Au collisions at ^JsNN — 200 
GeV have been recorded by the PHENIX experiment 
during Run4 (2003-2004). The much higher statistical 
precision of these data should allow an unambiguous re- 
sult on the important issue of charm flow. 
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